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In systems where the design inlet and outlet 
pressures P am b are maintained above the thermo- 
dynamic critical pressure P c , it is often 
assumed that heat and mass transfer are governed 
by single-phase relations and that two-phase 
flows cannot occur. This simple rule of thumb 
is adequate in many low-power designs but is 
inadequate for high-performance turbomachines, 
boilers, and other systems where two-phase 
regions can exist even though P^ m b > Pr. Heat 
and mass transfer and rotordynamic-f 1 uid-mechani c 
restoring forces depend on momentum differences; 
and those for a two-phase zone can differ signi- 
ficantly from those for a single-phase zone. By 
using a laminar, variable-property bearing code 
and a rotating boiler code, pressure and temper- 
ature surfaces were determined that illustrate 
nesting of a two-phase region within a super- 
critical pressure region. The method of corre- 
sponding states is applied to bearings with 
reasonaole rapport (within 5 percent for an 
oxygen/nitrogen example). 


1. INTRODUCTION 

As a rule of thumb the extent and nature of 
the cavity formed downstream of a conventional 
inlet and the minimum clearance in a bearing or 
seal can be controlled, even eliminated, by 
increasing the ambient pressure. Often systems 
are operated above the thermodynamic critical 
pressure in an effort to circumvent two-phase 
flow and cavity formation. Uncertainties in 
stability, load capacity, heat transfer, fluid 
mechanics, and thermophysical property variations 
are thereby avoided. 

Flow cavities are engendered by local depar- 
tures from classical thermodynamics, termed 
"metastabi 1 ities. " Classic among these are the 
work of Skripov [1] and applied work such as 
fracture, boiling, and condensing [2,3]. In many 
cases a two-phase flow cavity enhances bearing 
stability [4,5] and seal load capacity [6]. The 
inception of cavity flow is important to both 
turbomachine stability and load capacity. 

In bearings or eccentric shaft seals the 
fluid is compressed in the convergent geometry 
(maximum pressure is- reached before the minimum 
clearance) and expanded to subambient pressure. 


Because this pressure is often less than satura- 
tion, dissolved gases are released and flash 
evaporation can occur. Pressure recovery and 
recompression above ambient then follow. A 
simplistic version of the cycle is illustrated 
in Figure 1, for ambient pressures above the 
thermodynamic critical pressure. 

It is unusual to have a subcritical region 
nested within a supercri tical region without 
intervening physical barriers. Yet it has been 
demonstrated that two-phase flow cavities can be 
"nested" within a system with supercritical ambi- 
ent pressure and that conditions for the exist- 
ence of such nested cavities are prevalent for 
sharp-edge inlets and for eccentrically loaded 
high-performance turbomachines [7]. For the seal 
data of Figure 2 the working fluid was para- 
hydrogen (P c =1.2928 MPa, T c = 32.976 K, 
p c = 31.4 kg/irr) and the nonrotating shaft was 
fully eccentric. The pressure at the inlet to 
the third step was greater than P c . However, 
immediately downstream of the step the pressure 
on one "side" of the seal (e = 0) remained above 
P c while the pressure on the opposite "side" 

( e = it ) was nearly equal to the saturation pres- 
sure P s . Thus the ci rcumferenti al pressure 
profile encompasses all states from P-j n > P c 
at 9 = 0 to P $ at e = tt and represents a 
two-phase flow cavity nested within a super- 
critical pressure regime. 
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Figure 2 Pressure profile for fully eccentric, nonrotating, three-step 
space shuttle main engine (S SME) seal configuration. 

An embedded two-phase flow cavity zone was 
calculated for a highly eccentric bearing 
(Table 1) by using the constant-property Reynolds 
equations (no energy solution) for fluid hydrogen 
[ 7 J . There are no corroborating data. 

In this study we extend the investigation of 
nested subcritical pressure regions by using 
variable-property laminar conservation equations 
and apply the method of corresponding states to 
bearings . 


2. ANALYSIS 


The governing equations have been developed 
elsewhere [8] from the work of Dowson [9] and 
Fowles [10]. Only the basic form and assumptions 
are repeated herein. 

2.1 Assumptions 

1. The clearance/radi us ratio is so small 
that the cylindrical system can be unwrapped and 
represented by a system of Cartesian coordinates 
with smooth boundaries. 

2. The fluid is Newtonian, the flow is lam- 
inar, and inertia and body forces are neglected. 

3. Velocity derivatives du/dy and dv/dy are 
large compared with other gradients. 

4. Fluid properties and clearance are func- 
tions of coordinates (x,y,z). 

The most questionable assumptions are that 
the walls are smooth and that high shear flows in 
narrow passages can be described in terms of 
laminar flows well beyond the usual range of 


Reynolds numbers. In the authors' opinion, 
neglecting asperities and turbulence results in 
underprediction of viscous effects and would 
exacerbate the cavity zone formation. The 
consequences remain to be explored. 

2.2 Governing Equations 


The generalized Reynolds equation in the 
unwrapped Cartesian coordinates (x = circum- 
ferential, z = axial) becomes 
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Tne energy equation becomes 


(i) 
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where e = cT and cp is viscous dissipation. 

2.3 Solution Method 

The Reynolds equation was solved first for 
the pressure and velocity fields. The properties 
were evaluated at each grid point as a function 
of pressure and temperature. With the new set of 
properties the energy equation was solved. The 
process was iterated until the pressure and tem- 
perature convergence criteria were met. The 
grid used was 11 axial by 33 circumferenti al . 
Doubling the grid did not have a significant 
effect on the computations, and continuity was 
checked at each point. 


3. RESULTS 
3.1 Bearing 

The axial pressure gradients for a bearing 
are usually smaller than those for a seal, but 
the ci rcumferenti al gradients are usually much 
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larger. Both seals and bearings can be designed 
to take advantage of these pressure gradients to 
provide significant damping and stiffness, which 
enhance stability [11]. 

The pressure and temperature surfaces of 
Figure 3(a) represent the constant-property solu- 
tion of the conservation equations for a rota- 
tional speed of 90 000 rpm for the bearing 
configuration of Table 2. The cavity pressure 
was significantly lower than the saturation pres- 
sure although the ambient pressure was above the 
critical pressure. A similar result was found, 
without heat transfer, for an ambient pressure 
about eight times the critical pressure [7], and 
the bearing was highly eccentric. 

As the fluid approaches the minimum clear- 
ance, the temperature contours (Fig. 3(a)) become 
uniform, reflecting an increase in heat generation 
and axial convection out of the bearing. Immedi- 
ately downstream of the minimum clearance, bulk 
fluid is convected into the bearing to provide 
cooling but does not penetrate into the cavity 
zone. The heat generated becomes trapped, and 
the temperature continues to increase. 

Although solutions of the vari able-property 
equations (Fig. 3(b)) indicated only about an 
8-percent reduction in peak pressure and temper- 
ature, the two-phase zone was all but eliminated 
and was represented by a single point. 

Increasing the speed of rotation enlarged 
the two-phase zone. The pressure, temperature, 
viscosity, and heat capacity surfaces for this 
case are characterized in Figure 4. The property 
variations were significant, and the profiles 


were not symmetric because of the 0.1-MPa differ- 
ence in inlet and exit bearing pressures (Table 2). 

It is apparent that constant-property solu- 
tions tend to overpredict the load-carrying capac- 
ity and the extent of cavity formation. Note 
also that the thermal surface lagged the pressure 
surface and the energy was transported down- 
stream. This is not always the case; for ex- 
ample, when asperity contacts occur, the classic 
Blok problem enters and heat is transported up- 
stream (for fluid oil and low rotational speed) [12]. 

3.2 Extended Correspond! ng States 

Often for fluid (a) for a bearing or seal an 
excellent equation of state and data exist but 
not for fluid (b). In such cases it would be 
valuable to have a method of evaluating the seal 
or bearing performance with fluid (b). For the 
infinite journal bearing the theoretical load W 
supported by a stable constant-property 360° 
fluid film depends on the bearing parameter a 
[1,13] 
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Here the classic form is used rather than (1 - a^) , 
which provides more realistic results [14], and A 
is the surface area, A = 2irRL. 

Using the extended principle of correspond- 
ing states [15,16], it appears possible to pre- 
dict bearing behavior for single-component fluids 
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Figure 4 Analytic surfaces for a bearing (Table 2) rotating at 105 000 rpm. 


(and perhaps for multicomponent fluids) based on 
the known behavior of a reference fluid. Viscos- 
ity becomes 



Fluid 

V 

K 

P , 
c’ 

MPa 

Pc 3 

g/cnr 

M, 

g/g mole 

f 

h i 

F 

n 

n 2 

°2 

126.3 

154.78 

3.417 

5.082 

0.3105 

.4325 

28.016 

31.9988 

0.816 

1.22 

0.74 


and from equation (3) 

z x P 0 = Ph,/f = 1.5 x 1.22/0.816 = 2.242 MPa 

( ^1 T i 

A a ( p ,T) = Aq IP — , jj (5) Tq = T /f = 100/0.816 = 122.5 K 


where for one-component fluids 



Example : Let the reference fluid (subscript 0) 

be oxygen and the working fluid (subscript a) be 
nitrogen. The bearing temperature and pressure 
are 100 K and 1.5 MPa. Assuming constant proper- 
ties and 9 = $ = 1, solve equation (4) for 
nitrogen by using the thermophysical properties 
of oxygen. 

First, we need the thermodynamic properties 
at the critical point [17]. 


For comparison, the constant-property 
Reynolds equation was first solved by using the 
bearing parameter for nitrogen at P = 1.5 MPa 
and T = 100 K. Then the solution was repeated 
for the bearing parameter predicted for nitrogen 
based on oxygen properties at Pq = 2.242 MPa 
and Tq = 122.5 K. The reduced circumferenti al 
pressure profiles at the center of the bearing 
are illustrated in Figure 5 for nitrogen and for 
nitrogen predicted from oxygen properties. The 
agreement, within 5 percent, is better than 
expected as pressure predictions are less accu- 
rate than flow (continuity) predictions. This 
problem requires further exploration as the 
Reynolds equation is a form of the continuity 
equation. 

3.3 Boilers and Rotors 

A corollary class of nested two-phase flows 
occurs in the rotating boilers and cooled rotor 
windings of turbomachine systems. A rotating 
pressure containment vessel with axial flow 
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Figure 5 Reduced midbearing circumferential pressure profiles 
based on corresponding states for fluids nitrogen and oxygen. 


produces centrifugal and Coriolis body forces. 

For a rotating cylindrical boiler the radial 
position of the critical pressure surface can be 
adjusted to occur very close to the wall or the 
fluid-vapor interface; the two-phase region is 
nested. A generic boiler configuration is qiven 
in Table 3. 

Calculations for the generic boiler were 
carried out by using a rotary boiler code [18], 
with Freon-11 (P- = 4.41 MPa; T c = 471 K; 
p c = 0.553 g/cm 3 ) [19] as the working fluid. 

For the supercri tical pressure regime the heat 
flux follows the metastable Leidenfrost and film 
boiling loci, and the body force dependency fol- 
low the (0.25 to 0.35) power-law dependence [20]. 
Thus the prescribed wall heat flux becomes 

las _ n [0 - 25+0 - 02 i°9i 0 ("n (6) 

q ig 


The solution converged for wall pressures less 
than P c , but failed to converge for pressures 
greater than P c (conditions of Table 3). The 
problem was traced to the definition of Freon 
enthalpies in this region and requires a new 
equation of state. The properties code [17] is 
being integrated into the boiler code [18]. 

To the dynamicist the location and extent 
of the nested regions may become important only 
after the system instabilities are encountered. 


For the bearing considered herein constant- 
property solutions with no asperity contact tend 
to overpredict the load-carrying capacity and 
the extent of cavity formation, with the thermal 
surface lagging the pressure surface. Variable 
properties tend to smooth the surfaces and in 
some cases eliminate the two-phase cavity. In 
the authors' opinion, inclusion of turbulence 
would exacerbate the cavity formation, perhaps 
again approaching the constant-property case. 

The method of extended corresponding states 
appears to be a reasonable approach to predicting 
the behavior of a bearing lubricated with fluid 
(b) from the data and equation of state for the 
bearing with reference fluid (a). 

5. NOMENCLATURE 
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9c . 

h , h ( x , z ) 

L 

M 

n 

P 

Pr 


q 

R 

S 

T 

t 

Xi 

x 


y 

z 


n 


e 

x 

p 

n 


eccentricity ratio, e/C 
clearance, m 
specific heat, kJ/kg K 
diameter, m 
internal energy, kJ/kg 
eccentricity, m 

functions scaling working fluid to 
reference fluid 

acceleration (body force), m/s 2 

standard acceleration (body force), m/s 2 

dimensionless film thickness, y/C 

seal or bearing length, m 

molecular weight 

acceleration ratio, g/g c 

pressure, MPa 

reduced pressure, P/P c 

heat flux, W/nr 

radius, m 

entropy, kJ/kg K 

temperature, K 

time, s 

coordinate (x/L, z/L) 
circumferential coordinate, m 
clearance coordinate, m 
axial coordinate, m 
viscosity, Pa s 
shape functions 

circumferential coordinate, rad 
thermal conductivity, W/m K 
density, kg/m 3 
angular velocity, rad/s 


Subscripts 


c thermodynamic critical point (also 

superscript ) 
s saturation 

0 reference fluid 

a working fluid 
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transport within the near-critical region 
embedded within the fluid layer is largely 
unexplored. 
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TABLE 1. ANALYTICAL BEARING PARAMETERS 


Geometric 

Diameter, mm (in.) . . . 25.4 (1) 

Length, mm (in.) 25.4 (1) 

Clearance,, mm (in.) .... 0.0254 (0.001) 

Eccentricity 0.97 

Hydrodynamic 

Inlet pressure, MPa (psia) . . 13.8 (2000) 
Exit pressure, MPa, (psia) . . 13.45 (1950) 
Maximum pressure, MPa (psia) . 27.16 (3938) 
Minimum pressure, MPa (psia) 0.079 (11.5) 
Inlet temperature, K (°R) .... 20 (36) 

Speed, rad/s (rpm) 7850 (75 000) 


TABLE 2. BEARING PARAMETERS 


Material Steel 

Fluid Parahydrogen 

Geometric 

Length, mm (in. ) 76 (3) 

Housing thickness, mm (in.) . . 12.5 (0.5) 

Shaft radius, mm (in.) 37.5 (1.5) 

Radial clearance, mm (in.) .0.0254 (0.001) 

Eccentricity 0.7 

Hydrodynamic 

Inlet pressure, MPa (psia) . . . 1.72 (250) 
Exit pressure, Mpa (psia) . . . 1.62 (235) 
Inlet temperature, K (°R) . . 22.1 (39.8) 

Housing temperature, K ( R) . . 55.6 (100) 
Heat transfer coefficient, 

W/cm 2 K (Btu/hr ft 2 F) 

Shaft to fluid . 0.057 (100) 

Housing to fluid 0.0028 (5) 

Speed, rad/s (rpm) 9425 (90 000) 


10 996 (105 000) 


TABLE 3. ROTATING BOILER PARAMETERS 


External fins, mm (in.) 13 (0.5) 

Wall thickness, mm (in.) 13 (0.5) 

Container radius, mm (in.) 100 (4) 

Vapor interface radius, mm (in.) . . . 25 ( 1 ) 

Length, mm (in.) 300 (12) 

Speed, rad/s (rpm) 576 (5500) 

Wall pressure, MPa (psia) .... 4.62 (670) 
Inlet pressure, MPa (psia) .... 3.1 (450) 
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